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Abstract 

It is well-known that a delta potential well in ID has only one bound state but that in 3D it 
supports an infinite number of bound states with infinite binding energy for the lowest level. We 
show how this also holds for the less familiar 2D case, and then discuss why this makes 3D delta 
potential wells unphysical as models of interparticle interactions for condensed-matter many-body 
systems. However, both 2D and 3D delta wells can be "regularized" to support a single bound level 
which in turn renders them conveniently simple single-parameter interactions, e.g., for modeling the 
pair-forming dynamics of quasi-2D superconductors such as the cuprates, or in 3D of other supercon- 
ductors and of neutral-fermion superfluids such as ultra-cold trapped Fermi gases. 

PACS numbers: 03.75.Ss; 03.65.-w; 03.65. Ge; 74.78.-w 
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Resumen 

Es bien sabido que un pozo de potencial delta en ID tiene un solo estado ligado pero que en 3D 
tiene un niimero infinito de estos estados con una energia de "amarre" infinita para el nivel mas bajo. 
Aquf mostramos como esto tambien ocurre para el caso bidimensional que es menos familiar, para luego 
discutir por que los pozos de potencial delta en 3D no son ffsicos como modelos de interacciones entre 
particular para sistemas de muchos cuerpos en materia condensada. No obstante, ambos pozos delta en 
2D y 3D pueden ser regularizados para soportar un solo nivel ligado lo cual los convierte convenientemcntc 
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en interaccioncs dc un solo parametro, por ejemplo, para modelar la dinamica de formation de pares en 
supcrconductores casi-bidimcnsionales tales como los cupratos, o en 3D la formacion de pares en otros 
superconductores y en superfluidos fcrmionicos neutros tales como los gases de Fermi atrapados ultrafrios. 
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I. INTRODUCTION 

The study of physical systems in dimensions lower than three has recently shed its purely academic 
character and become a real necessity to describe the properties of novel systems such as nanotubes [T] , 
quantum wells, wires and dots 2, 3], the Luttinger liquid PUS], etc. Reduced dimensionality describes 
superconducting phenomena in quasi-2D cuprates where pairing between electrons (or holes) is essential 
[5]- Whatever the actual interaction between two electrons (or holes) in a cuprate might ultimately turn 
out to be, the attractive delta potential is a conveniently simple model to visualize and to account for 
pairing, an indispensable element for superconductivity and neutral-fermion superfluidity. It enormously 
simplifies calculations. Bound states in a delta potential well in ID and 3D are usually discussed in 
textbooks, but not in 2D. Refs. (i and [7] discuss this from a more rigorous mathematical viewpoint 
without explicitly solving the Schrodinger equation, e.g., for the bound energy levels. Here, this gap 
is filled by analyzing the 2D time-independent Schrodinger equation with a delta potential well that is 
then "regularized" [5] to reduce its infinite bound levels to only one. The single-bound-level case suffices 
for now since, e.g., the well-known simple Cooper/BCS model interaction [5] mimicking the attractive 
electron-phonon pair-forming mechanism, but requiring two parameters (a strength and a cutoff) instead 
of the regularized ^-potential's only one (a strength), can also be shown to support a single bound state 
|10j in the vacuum or two-body limit. Were it not for the (momentum-space) cutoff parameter, the 
Cooper/BCS interaction in coordinate space would also be a (5-potential, and indeed becomes such as the 
cutoff is properly taken to infinity. 

From elementary quantum mechanics we first recall the bound-state energies E < in a potential 
"square" well of depth Vo and range a, a common textbook example studied in ID [11] and 3D [12]. In 
ID the ground-state energy Eq of a particle of mass m can be expanded for small V$a as 



Thus, in ID there is always at least one bound state no matter how shallow and/or short-ranged the well. 
Similarly, in 3D for a spherical well, an expansion of Eq in powers of r\ = Vqcl 2 — h 2 ir 2 /8m > gives 



Thus, in contrast to ID, a minimum critical or threshold value for V$a 2 of fi 2 7r 2 /8m is needed in 3D 
for the first bound state to appear. Clearly, both ID and 3D cases are perturbative expansions in an 
appropriate "smallness" parameter, V$a or 77. As in ID, a 2D circularly symmetric well of depth Vo and 
radius a always supports a bound state, no matter how shallow and/or short-ranged the well. However, 
this instance is non-perturbative as it gives |13] for the lowest bound-state energy 
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which cannot be expanded in powers of small Vqo 2 
essential singularity at A = 0. 



since it is of the form /(A) 



= e 



l/A 



0, i.e., has an 
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In this paper we discuss how, just as in the better known 3D case, the 2D potential well — voS (r) , 
vo > 0, also supports an infinite number of bound states with the lowest bound level being infinitely 
bound for any fixed vq. For an N — * oo many-fermion system interacting pairwise via a delta potential, 
arguments based on the Rayleigh-Ritz variational principle show that the entire system in 3D would 
collapse to infinite binding energy per particle E/N —> — oo and infinite number density n = N/V — > oo. 
This occurs since the lowest two-particle bound level in each J-well between pairs is infinitely bound, 
for any fixed vq. To avoid this unphysical collapse one generally imagines square wells in 3D (and also 
in 2D) "regularized" [8] into 6 wells —vqS (r) that support a single bound-state, a procedure leaving an 
infinitesimally small vq. The remaining ^-potential well is particularly useful in condensed-matter theories, 
e.g., of superconductivity [14] or neutral-fermion superfluidity [151116) . where the required Cooper pairing 
can arise |17j from an arbitrarily weak attractive interaction between the particles (or holes) . 

After beginning with a d-dimensional expression for the delta potential in Sec. II, we summarize how 
bound states emerge in ID and 3D by recalling textbook results. In Sec. Ill we analyze in greater detail 
the less common 2D problem. In Sec. IV we sketch the use of "regularized" 2D S potential wells for 
electron (or hole) pairing in quasi-2D cuprates and in Sec. V we give details for the 2D case. Sec. VI 
offers conclusions. 



II. REVIEW OF DELTA POTENTIAL WELLS IN ID AND 3D 

The attractive square potential well in d dimensions 

V(r) = -V 9(a-r), (4) 

where the Heaviside step function 9{x) = i[l + sgn(a;)], a is the well range, and Vo > its depth. An 
attractive delta potential — v S(r) (v > 0) can then be constructed from the double limit 

-V d(a-r) > -v Q S(r), (5) 

Vo — > oo, a — > 

3 a d V = const. 

where vq = Cda d Vo is a positive constant, with Cd = ir d / 2 /T (d/2 + 1) as follows on integrating both sides 
of (O over the entire d-dimensional "volume" and recalling that J d d r5(r) = 1. We seek the bound-state 
eigenenergies E < from the time-independent Schrodinger equation for a particle of mass m in potential 
dU), namely 

V 2 *(r)-^[V(r) + |^|]*(r) = 0, (6) 
nr 

where E = — \E\ . 

In ID the solutions of J6]) (with r > taken as |x|) for x ^ are ^>(x) — e ±px . These functions have 
a discontinuous derivative at x — a in the delta potential limit ([5]) where 

lim 2aVo = vo < oo, (7) 

Vo — ^oo, a — >0 
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and there is always a (single) bound-state energy E — —mv 2 l /2h 2 (vq ^ 0) [TT]. Note that in the integral 
method of Ref. T3] applicable to shallow wells where \E\ <C max |V(x)| , E for ID would be given as 

i 2 



E 



m 
2h? 



dxV(x) 



which for V(x) 



-VqS(x) becomes 
E = 



2ft 2 



dxS(x) 



TTIVq 

lift 2 "' 



(8) 



(9) 



which agrees with Ref. [TT] and is consistent with (jTJ) . The ID <5-potential well has proved very convenient 
in modeling [181 119] self-bound many-fermion systems in ID, and in understanding Cooper pairing [20) 
as well as the BCS theory of superconductivity [21] . 

For the potential ^ in 3D the particle wave function in spherical coordinates [22j is ^(r) = Ri(r)Y[ m (9, < 
where (Ref. [5^, p. 722) Yi m (9,4>) are the spherical harmonics and Ri{r) the radial wavefunctions. For 
< r < a the finite (or regular) radial solutions are spherical Bessel functions of the first kind ji(Kr) 
of order I, with K 2 = 2m(Vb — \E\)/h 2 , since ji(Kr) < oo at r = 0. For r > a the linearly-independent 
radial solutions are the so-called modified spherical Bessel functions ki(kr), with k 2 = 2m \E\ /h 2 , where 
ki{kr) decays exponentially as r — > oo. The boundary conditions at r — a expressing the continuity of 
the radial wave function Ri (r) and of its first derivative can be combined into the single relation 



dji (Kr) I dr 



Ji(Kr) 



dki(kr)/dr 



ki(kr) 



(10) 



Taking I = and recalling (Ref. |23|, pp. 730, 733) that jo(x) — sin(x)/x, and ko(x) = e x /x, (fT0|) gives 

Kcot(Ka) = -k, (1 = 0). (11) 

We can write the I = bound-state energies E n = —h 2 /2ma 2 (7r 2 /4 + e 2 ) , where e n are the dimensionless 
roots of (fTTj) . with n = 1,2, .... The standard graphical solution [T^] of condition (jTTJ) shows that there 
are precisely n bound I = states whenever the well parameters are such that [24] 

'2mV a 2 \ 1/2 



(n-l/2)7r < 



h 2 



< (n + l/2)7r; (n = l,2, 



(12) 



Thus, the first bound state (n — 1) appears when Vqcl 2 > n 2 h 2 /8m, as was mentioned below ([2]), and 
n = 2 requires a deeper well depth Vb and/or larger well range, etc. 

The 3D delta potential well — vq<5 (r) , as defined in ([5]), integrated over all space gives 



vq 



d tvoS(t)= lim 

Vb — >oo, a — >0 



d 3 rV 9 (a - r) 



4tt.. 3 
hm — V a 6 



< 



(13) 



Hence, as Vb — > oo, a — > the middle term in p^|) (2mVba 2 /fi 2 ) 1 / 2 = (3mvo/27rh 2 a) 1 ^ 2 — ► oo, so that the 
number of bound-states n in the 3D delta potential well — vq8(y) is infinite for any finite fixed strength vq. 
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III. 2D DELTA POTENTIAL WELL 

This same result holds in 2D but is not as apparent. Here the solutions of © are $(r) = f(r)e w ^, 
with v = 0,1,2,... and the angular variable — 7r < <j> < tt. For < r < a the radial solutions which 



are finite at r — are cylindrical Bessel functions J v (Kr) = \J 2Kr / nj„_ 1 / 2 (Kr) (Ref. [23], p. 669) of 
integer order is, with K 2 = 2m(Vo — | £7 1 ) / 7i 2 . For r > a, as linearly-independent solutions one has the 
modified Bessel functions K v {kr) with k 2 = 2m \ E\ /ti 2 , which are regular as r — > oo. The two boundary 
conditions at r — a can again be written as a single relation 



dJ u (Kr)/dr 



dK v {kr)/dr 



(14) 



As we want to ensure against collapse in our many-body system interacting pairwise with the 6 poten- 
tial, it is enough to show this for the lowest bound level with v = 0. In this case (jT4"]) becomes, since 
dJ (Kr) /Ar = -KJ X (Kr) and dK (kr) /Ar = -kK x (kr) (Ref. [25], p. 361 and 376, respectively), 

Ka^± ka^± (15) 
Ka Jo(Ka)- ka K (kaY (15) 

Since K\(x) > Kq(x) > for all x, the rhs of (|15p is always a positive and increasing function of fca; 
it is plotted in Fig. 1 for Vq/ \E\ = 300 (dashed curve). As for the lhs, Jo (x) oscillates for all x so 
that it diverges positively whenever Jq(x) = 0, then changes sign and thus drives the lhs to — oo (see 
full curve in figure). Clearly, there is always an intersection (bound state, marked by dots in figure) 
between two consecutive zeros of Jo (a;). For a given interval in ka, the closer these poles are, the more 
bound-states there will be. Thus, for any given square well, all of the allowed bound-states lie inside an 
interval between and k max a, where fc max = (2mVo/h 2 ) 1 ^ 2 . In such an interval the number n of bound 

1/2 

states (zeros) will be n = INT(aA: max a/7r), with a = (Vq/ \E\ — 1) ' where the INT(x) function rounds a 
number x down to the nearest integer. Of course, the expression for n is only valid after the appearance 
of the first pole. Then for Vq/ \E\ = 300 as in Fig. 1, n = 3 in the interval between and k max a — 0.5. 
In Fig. 2 are shown the bounds for Vq/E = 2700 where there are n = 7 bound states in the interval 
between and fc max a = 0.4 as it should be. 

To construct a delta potential well —vo6(r) in 2D from the finite-ranged well (|4]), and through ([5]) 
ensure that J d 2 r5 (r) = 1, requires that 

lim Voira 2 = vo < oo. (16) 

Vo — >oc, a — >0 

Thus, as long as \E\ is finite ka = (\j2m \E\ /h 2 )a — > and we can use (Ref. [17] . p. 612) 

Vo — >oo, a — >0 

xKi(x)/ Kq(x) — ► — 1/lnx for the rhs of (fl~5|) . In this case the number n of bound states for 8- well 
corresponds again to the number of zeros of the lhs of (fl~5f but in the delta limit. Here, from (fT6|) 
Ka = ka (Vq/ \E\ — l) 1 ^ 2 ► J 2mvo I < oo (not necessarily 1). We will see below that the 

Vo — >oo.a — >0 

case Ka <C 1 corresponds to the shallow 2D potential well of Ref. [13]. But even if Ka is not <C 1, Bessel 
functions oscillate for large argument although their period is not constant. In this latter case (Ref. [2"5] . 
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p. 364) J\(x) — > \j2fwx cos(a; — 37r/4) allows locating the zeros of the lhs of (fT"5|) which as VqI \E\ is 

a:»l 

increased approach each other on the ka axis, so that in the delta well limit as Vq/ \E\ — ► 00 the number 
n of bound-states increases indefinitely. Moreover, rewriting (|15[) as 

aSnJ^ainJ^o^n) - a^Ai (x„ ) J Q (ax n ) = 0, (17) 



bound states are easily identified from Fig. 3, where the roots of (fT7)) . say x„ = k n a — ( -y/Im /?i 2 
are seen to form an infinite set as a — > 00. Therefore the 2D delta potential well supports an infinite 
number of states, for any fixed Vq, precisely as in the 3D case, this being the main conclusion of the paper. 
Table 1 shows the first few (numerical) x„ roots where E n = —h 2 x 2 l /2ma 2 , for three extreme values of 
Vol \E\ . 

Applying the integral method of Ref. [T3] for v = for a shallow potential well, i.e., Vo — > and 
\E\ <C Vo, one can take both Ka and ka — > 0. Thus, we can use xKi{x)/Kq(x) — ► — 1/lnx in the rhs 
of (HID, and in the lhs of dUD we note that (Ref. [25], p. 360) J v {x) — > x v /2 v v\, with v = and 1, so 
that xJi(x) I Jq(x) — ► x 2 /2 . Hence we write (fT5j) as 

1 lnfca 

"(i^^^"' (18) 
so that on putting Vq — \E\ w Vq (| 18[) becomes precisely <J3J> - In fact, for any shallow 2D circularly- 
symmetric potential well V (r) , the first bound state in Ref. [13] is given by 



2 



E ~ ^ exp ( —H z /m 

2ma l 



drrV(r) 



(19) 



which for potential |4} reduces to ([3]). This result in the delta limit of (fT6| finally becomes 



ft 2 

E ~ ex P (-2» 2 7r/m«o) , (20) 



where w < 00. 



IV. NEED TO REGULARIZE IN CONDENSED-MATTER SYSTEMS 

Real condensed matter systems are made of many particles (bosons and/or fermions) interacting via 
attractive and/or repulsive forces. Attractive forces between fermions can form pairs needed for many 
properties such as superconductivity in solids or superfluidity in fermion liquids or trapped atomic fermion 
gases. However, addressing these problems with a physically realistic interaction is oftentimes difficult. 
As in ID with a "bare" 5-potential well, a regularized attractive <5-well prevents collapse in 3D, provides 
the required pairs in either 2D or 3D and, of course, simplifies calculations. 

It is easy to imagine a trial wave function whereby, with an attractive bare 5-function interfermionic 
interaction (i.e., before regularization), a 3D N- fermion system would have infinitely negative energy- 
per-particle (as well as infinite number-density). This is because the lowest bound level of the two-body 
(5-well is infinitely deep in 3D, and indeed also in 2D, as was shown in the preceding sections. By the 
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Rayleigh-Ritz variational principle the expectation energy associated with the trial wave function is a 
rigorous upper bound to the exact iV-fermion ground-state energy, and hence produces collapse of the 
true 3D ground state of the system as N — > oo. In this picture, each particle "makes its own well" but 
will attract to itself every other particle, two for each level, to minimize the trial expectation energy. We 
thus get an iV-fermion system as schematically sketched in Fig. 4 (where the Pauli exclusion principle 
is explicitly being applied) that collapses as N — > oo. To avoid this unphysical collapse in 3D, and at 
the same time ensure pair formation in either 2D or 3D, one can "regularize" [5] the 2D and 3D finite 
interparticle potential wells so that in the limit the corresponding <5-well possesses only one (s-wave) 
bound state. This also occurs with the Cooper/BCS model interaction [TH] definable in any d and with 
the bare (5-well in ID. The single-bound-state (5-well then ensures that only pair "clusters" form, in agree- 
ment with quantized magnetic flux experiments in either elemental (2f5j 27J or cuprate superconductors 
[28] in rings where the smallest flux trapped is found to be h/2e (with h being Planck's constant and 
e the electron charge). This contrasts with h/e as London conjectured just on dimensional grounds, as 
well as with h/ne (n = 3, 4, ...) which is not observed in superconductors, as one would expect in vacuo 
in other many-particle systems with attractions that produce clusters of any size. The fact that only 
pair clusters occur with electrons (or "holes" ) in superconductors is likely associated with clusters forming. 



V. REGULARIZED 2D DELTA WELL 

Regularization in either 2D or 3D starts from a finite-range square well and yields (5-well with an 
infmitesimaHy small strength vq, as we now illustrate. To be specific, we concentrate on the regularization 
of the 2D finite potential wells needed to mimic, in a simple way, the presence of Cooper pairs in 
superconductors. We thus seek a two-body square well interaction such that, in the <5 limit, it possesses 
only one (s-wave) bound state. Following Ref. f8J, for d — 2 we substitute ||3J| by an effective two-body 
square- well interaction V a {r) which in the limit a — > + becomes — voS(r) with Vq > 0, and is given by 

V a (r)= H \ 2 M a-r), (a > a > 0) (21) 

2m a z in |a/ao| 

where m* = m/2 is the reduced mass of the pair, a is still the well range, and eto is an arbitrary parameter 
that measures the actual strength of the interaction. Potential (f2"Tj) in the delta limit © gives a 5-well 
strength vo(a) = —jd rVa(r) = — irh 2 /m* In | a/ao ] > 0, and is thus infmitesimaHy small as a — > + . 
However, this parameter can be eliminated in favor of the binding energy B2 > of the single level, which 
now serves as coupling parameter. Indeed, using v (a) in equation (|20[) for a shallow-well as in Ref. [13] 
but with to* instead of m, we obtain for the lowest energy 

E = —^ = -B 2 , (22) 

7710q 

where < B2 < 00 is the magnitude of the pair binding energy. This straightforward procedure then 
guarantees a simple finite-lowest- energy level well — as in ID, see (9). Once we set the regularized two- 
body interaction model, result (|22[) can be varied as the coupling describing our superconductor model. 



One possibility we have now is to fix B2, which fixes the value of ao; the second possibility is to set B2 
by fixing the parameter ao, which can represent the range of the wave function of the particles. One can 
also reduce the infinite number of bound states to only one by shifting the center of the 2D (5-potential 
from the origin along the radial axis |29] ; however the topology of this new ^-potential is unsuitable to 
simulate real interactions between electrons. 

In 3D regularization proceeds similarly [5] as in 2D except that there is no log term in (|2"Tj) , and 
instead of the binding energy (|22[) as coupling parameter one employs the s-wave scattering length which 
is well-defined even if, unlike 2D, the well is too shallow to support a bound state. 

In an A-fermion system interacting pairwise via a regularized 5-potential, fermions in the Fermi sea 
bind each other by pairs only, as required by magnetic flux quantization experiments [26 27|[28], see 
Fig. 5. The (5-potential has been used extensively in the literature [301 131] to mimic the pair- forming 
interfermion interaction in, e.g., quasi-2D cuprate as well as in otherwise 3D superconductors [31U32j and 
neutral- fcrmion superfluids [151 116] . 

VI. CONCLUSIONS 

A graphical proof was provided of how a 2D (5-potential well supports an infinite number of bound- 
states as does the more familiar 3D (5-potential well. Using Rayleigh-Ritz variational-principle upper- 
bound arguments, we then illustrated how in 3D the binding energy-per-particle of an A-fermion system 
must grow indefinitely as the number of fermions increases. In order to prevent this unphysical collapse 
in modeling such a system one can use regularized (5-potentials in 3D as well as in 2D that by construction 
support a single bound state. This provides useful interfermion interaction models to study 2D and 3D 
condensed-matter problems. An appealing motivation for regularized (5-potentials is that they fit easily 
within the framework of the time-independent Schrodinger equation in either coordinate or momentum 
space. Indeed, solving the two-body problem in the Fermi sea allows one already to exhibit Cooper 
pairing phenomena which is the starting point for any treatment of superconductivity or neutral-fermion 
superfluidity. 
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Tables 



Vol \E\ 


X\ 


X2 


X3 


X4 


10 


0.3738 


1.4216 


2.4674 


3.5137 


10 3 


0.0218 


0.1248 


0.2245 


0.3240 


10 5 


0.0123 


0.0223 


0.0323 


0.0422 



Table 1: First few roots x n — (i/2m \E n \ /fr 2 )a of (TTT)) for bound-states E n of 2D potential well 
according to (|15p . for different values of Vq/ \E\ . 
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Figure 1: Rhs (dashed curve) and lhs (full curve) of (fT5|) for the 2D well with Vq/ \E\ = 300. Intersections 
of both curves marked by dots signal bound states. There is always a bound state between every two 
consecutive zeros of Jo(x), or poles of the lhs of (|15| . 
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Figure 2: Same as Fig. 1 but for Vq/ \E\ = 2700 suggesting that the number of bound states increases 
indefinitely as the potential well approaches the Dirac S- well limit ([5]). 
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Figure 3: The bound states for the 2D (5-well are associated with the zeros of the lhs of equation (|T7|) . 
plotted on the vertical axis. For Vq/ \E\ = 1000, this graph illustrates the roots (bound states) of (jTTJ) 
which become an infinite set as we approch to the (5-well limit. 
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2D or 3D delta 

potential 

(schematic) 




Figure 4: An iV-fermion system with the <5-well pairwise interactions produces collapse as N — > oo in 3D 
since both the binding energy per particle and the particle density diverge. 
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dimer gas 
(schematic) 




Figure 5: A dimer gas formed by single-bound-state regularized (5-wells, schematically depicted. In this 
case, the 3D many-fermion system will not collapse since the Pauli exclusion principle prevents more than 
two particles from being bound in a given well. 
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